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Local quenches in frustrated quantum spin chains: global vs. subsystem equilibration 
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We study the equilibration behavior following local quenches, using frustrated quantum spin chains as an 
example of interacting closed quantum systems. Specifically, we examine the statistics of the time series of 
the Loschmidt echo, the trace distance of the time-evolved local density matrix to its average state, and the 
local magnetization. Depending on the quench parameters, the equilibration statistics of these quantities show 
features of good or poor equilibration, indicated by Gaussian, exponential or bistable distribution functions. 
These universal functions provide valuable tools to characterize the various time-evolution responses and give 
insight into the plethora of equilibration phenomena in complex quantum systems. 



I. INTRODUCTION 



The equilibration properties of closed quantum systems 
have been a topic of growing interest. As opposed to equili- 
bration in open quantum systems coupled to an infinite bath, a 
closed finite quantum system does not admit asymptotic fixed 
points, unless the initial state is an eigenstate of the evolution 
Hamiltonian. Except for this trivial case, the evolution state 
p (t) oscillates without converging to any limiting point. On 
the other hand, equilibration in closed quantum systems may 
be defined by resorting to probabilistic notions |[l]-|gl- Specif- 
ically, one says that an observable equilibrates if it remains 
close to its average for most of the times along the evolution 
trajectory. Alternatively, equilibration takes place if the fluctu- 
ations around the average are small. For typical initial states, 
the authors of l|3l|4[] proved a number of bounds on the fluctua- 
tions that imply equilibration, provided the subsystem is small 
compared to the environment. Here "typical" means "taken 
uniformly at random from the space of pure states". The in- 
sightful results in fj, ^ leave nonetheless open the question 
of equilibration for a number of interesting cases. Namely 
when sizes are such that the environment cannot be consid- 
ered overwhelmingly larger than the subsystems and/or when 
the initial state cannot be considered typical. Both such fea- 
tures are present in the experimentally relevant situation of 
small quenches on a finite system liiSJ. This is the situation 
studied in this article. More precisely, we suddenly change the 
evolution Hamiltonian only on a subsystem S and then mon- 
itor equilibration features of the subsystem as well as of the 
total system by studying the full statistics of global (i.e. with 
support on the whole system) as well as of local observables. 
For the local subsystem we choose the subsystem magnetiza- 
tion as a physically natural observable. To obtain informa- 
tion on all possible observables defined on S we monitor the 
trace distance of the subsystem density matrix from its equi- 
librium value. As a global observable we take the Loschmidt 
echo (LE). The Loschmidt echo has the additional appealing 
features that its time averages provides a bound on the fluc- 
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tuations on any sufficiently small subsystem. However, in the 
situations under study most of the bounds are too loose and 
equilibration must be checked by an explicit computation. 

To be more specific, the idea is the following. We prepare 
the system at i = in the ground state of the Hamiltonian 
Ho. At t > we suddenly modify the Hamiltonian in the 
region S. The system then evolves without further disturbance 
with Hamiltonian H = Hq + Vs, where the term Vs only 
acts on subsystem S. The equilibration properties are then 
evaluated by monitoring local quantities of S as well as global 
quantities of the entire system. Throughout this paper we will 
use the terms "local" resp. "global" for quantities which refer 
to subsystem S resp. the entire system. The idea is to study the 
interplay between equilibration of these "local" and "global" 
quantities. 

After a sufficiently long evolution time, the trajectory p [t) 
defines an emerging equilibrium average ensemble given by 
'p :— p [t) := limT^oo T^"^ J^ p (t), where the observation 
time T has been sent to infinity. The average density matrix p, 
in the eigenbasis of the evolution Hamiltonian H, has the form 
of adephased state: p = J2nP'n\''^){M' where p„ = KtiIt/io)! 
(iV'o) initial state), and for simplicity we assume here that the 
eigenenergies En are non-degenerate flS- In the same fash- 
ion, if we monitor an observable O over a long time, its av- 
erage value will be given by O :~ {O (t)) :— tr [p (t) O]. In 
finite systems we can safely bring the time average inside the 
trace and obtain O — ti [pO] . The equilibration properties of 
an observable are encoded in its variance, or more precisely in 
its full statistics. The quantity {O (t)) is seen as a stochastic 
variable with the uniform measure dt/T over the observation 
time interval [0, T] (T will be sent to infinity). We can say that 
the observable O has good equilibration features if its proba- 
bility distribution Pq shows concentration around its mean 
value. Explicitly, the probability distribution Pq is given by 
Po (x) = 5 {{O (t)) — x). To fix ideas further, we can say 
that good equilibration coiTesponds to a small variance AO^, 

or more precisely to a large signal to noise ratio O/v AO^. 
Clearly Pq depends on the observable O but also on the initial 
state \4>q) and on the evolution Hamiltonian H. 



II. PRELIMINARIES 

The model we study here is the frustrated antiferromagnetic 
Heisenberg spin- 1/2 chain, with a magnetic field acting only 
on a subsystem S containing Ng adjacent sites. Its Hamilto- 
nian is given by 
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(1) 



where Ji and J2 denote the nearest-neighbor and next- 
nearest-neighbor couplings respectively, Ji is set to 1 here- 
after, Sj are spin- 1/2 operators, and periodic boundary con- 
ditions are applied throughout, i.e. Sj+n = Sj. The ini- 
tial, i.e. pre-quench, Hamiltonian Hq is given by Eq. ([1]) with 

hs = hg . The local field hg is then suddenly changed, such 
that the evolution Hamiltonian H is given by Eq. ([T]l with 



hs 



_ hif) 



The couplings Ji, J2 are held fixed during the 



quench process. 

If we consider the hard-core boson mapping of the spin 
variables given by S^ — h'.bj — 1/2, the local field hs is equiv- 
alent to a local chemical potential which confines the particles 
in region S. The Hamiltonian Eq. ^ commutes with the total 
magnetization S^^^^. Since in the numerical studies discussed 
below the initial field ^,5 is not chosen excessively large, we 
verified that the ground state of Hq (and so the evolved state 
too) always belongs to the sector with zero total magnetiza- 
tion S'^Qj = 0. This means that the system has N/2 hard core 
bosons which initially prefer to occupy the region S; still the 
hard-core constraint ((6J)^ = 0) does not permit to have more 
than one particle at the same site [23}]. 

The local field also breaks translational invariance. This en- 
larges the dimension dp of the ground state manifold, i.e. the 
vector space which contains the starting state, \4'q) and the 
evolved state \tjj {t)) = exp {—itH) \tj}o)- This feature is im- 
portant. In fact, do plays the role of an effective Hilbert space 
dimension. In systems with many symmetries, do can be sig- 
nificantly reduced with respect to the total Hilbert space di- 
mension. As a consequence, a system with large do has com- 
paratively smaller finite size effects compared to a system with 
the same number of sites but more symmetries and smaller do • 

To obtain information about the system's equilibration 
properties we consider various time dependent observables 
/ (i) and compute their statistics using the uniform measure 
dt/T over the interval [0, T], in other words we compute the 
probability distribution Pf of the random variable / {t) via 
Pf (x) := S (/ (t) — x). We are mainly interested in the con- 
centration properties of Pf, that is its "peakedness". 

To monitor what happens to subsystem S we first 
consider its magnetization rng {t) :— {Sg {t))/Ns — 
{ip [t) 1 5*11 V" (t))/Ns- As for any other observable O, it can 
be expressed in the basis of the eigenvectors of the evolution 
Hamiltonian in the following way: 

(0(t)> =0+J2 '2{n\0\m)c^CnCos[t {E^ ~ En)] . (2) 



Here c„ = (rilV'o)^ and we assumed that both the observable 



O and the initial state |i/)o) can be chosen real in the basis of 
H. 

The local magnetization Sg, although being quite natural, 
is just one of the 2^^' linearly independent observables which 
we can define on 5*. To obtain information on all possible ob- 
servables on S one has to consider the density matrix of the 
subsystem S: ps (t) — trg-lV-" {t)){'4' (0 I- A convenient quan- 
tity which encodes the equilibration properties of all observ- 
ables in S is then given by the trace distance between ^5 [t) 
and its equilibrium value pg, i.e. 

Dg{t):=\\pg{t)-J^\\^, (3) 

where the norm-1 of a matrix O is given by \\0\\-^ := 
^trVOUJ. The trace distance Dg in Eq. (O was first consid- 
ered in [4]. If the time average of Dg (t) is small, the evolved 
state of the subsystem pg (t) stays close to its average JJg for 
most of the time. This can be made quantitative by noting that 
Dg is a positive quantity and using Markov's inequality: 



Prob [Dg (t) >e]< 



D, 



valid for any positive e. This implies that, for any observable 
O defined on 5*, the time evolved quantity {O (t)) is close to 
its mean O for most of the time. That is if Dg is small, then 
one finds good equilibration for all observables in S. This can 
be shown by noting that, if the observable O has its spectrum 
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then H: 




\tr[pg{t)0] - tl[-p^O]\ < (o„,ax - 0„,in) Dg{t). (4) 

The authors of [MJ] also proved the following interesting in- 
equalities 



(5) 



Here dg is the dimension of subsystem S, equal to 2 ^' in our 
case, the effective dimension of a density matrix d°^ (f) is de- 
fined as the inverse purity, i.e. d°^ {a) = 1/tr (c^)- Finally, 
/>g is the equilibrium state restricted to the complementary of 
5, i.e. its environment S. Especially the second inequality 
of Eq. (|5]l is very useful. It gives a bound on Dg by resort- 
ing to a global quantity d°^ (p). As was shown in 141191 UOtl . 
for large quenches 1/d'^^ (p) is typically exponentially small 
in the system size: l/d^^ (p) « cxp (—const, x N). There 
is, however, an interplay between system size and quench 
amplitude. For moderate quenches, and/or when the sys- 
tem size is not too large, the constant appearing in the ex- 
ponential can be very small. This can be seen using per- 
turbation theory for small quenches which predicts i9l-lllll 
l/d°^ (p) « exp (—const, x (5/i| x Ng). Thus, when Ng is 
not large compared to (Shg)^^, the second bound in Eq. ^ 
is not effective, since Dg < 1 by definition. Does this mean 
that equilibration in this case will not take place in subsystem 
SI This is precisely the question we are going to address here. 
Note that in our simulations with Ng — 4 the second bound in 
Eq. (|5]) is not effective as long as 1/d'^^ (p) > 0.015, whereas 



in most of our simulations we have 1/(F^ [p) > 0.9 (see be- 
low). 

As a prototype quantity to monitor global equilibration of 
the entire system, we study the Loschmidt echo (LE) Ill2l - t20ll 

£(t) = |(V;o|expHiJf)|V'o)|'. 



Here H is the evolution Hamiltonian, and |V'o) is the initial 
state. C (t) measures the probability that the state \tp (t)) 
evolves back to the pre-quench initial state |^o)- In this sense, 
the LE is the expectation value of a particular operator, i.e. 
the projector onto the initial state \ipo) {ipo\- This operator has 
clearly support on the entire system, and therefore the LE is 
the time dependent expectation value of a "global" quantity. 
The LE has many interesting interpretations (see e.g. 191 UM 
and reference therein). Most important for us, the time av- 
eraged LE is precisely the inverse of the effective quantum 
dimension defined above. In fact, writing the LE in the eigen- 
basis of H one finds 
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Figure 1 : Lowest 5 energy levels of Hamiltonian Eq. (TJ on a ring of 
A^ = 16 sites for different coupling ratios J2/J1 as a function of the 
field hs applied on A''s = 4 contiguous sites. 



III. FIELD-ENERGY-DEPENDENCE AND QUENCHES IN 
DIFFERENT REGIMES 



£ (t) = £ + 2 V PnPm COS [t {En - E^)] 



(6) 



where £ - ZnPl - tr (p|V'o)(^o|) = tr [(p)^] = 
1/(P^ (p). Given its simple form in terms of the weights 
Pn — knl ' '^he probability distribution of the LE is closely 
related to the distribution of the p„. 

Once again, because of the bounds in Eqs. ^ and (|4]l, a 

sufficiently small average of the LE implies equilibration in 

all (sufficiently small) subsystems. Moreover, using the bound 

£" < nlC proved in [9]] with n = 2 we obtain that the 

2 

variance satisfies AC^ < C , and thus a small average implies 

small variance, with a signal to noise ratio greater or equal 

to 1. This scenario can be considered the standard road to 

equilibration. The converse need not be true, i.e. one can have 

a small variance for the LE but a large mean. The possibility 

of such alternative scenarios is what we are going to explore 

in this paper. 

Details of the numerical computation: In the numerical 
computations we take advantage of conservation of total mag- 
netization to reduce the Hilbert space dimension. To find the 
initial state |^o)> we first diagonalize Hq in the lowest mag- 
netization sectors, S^^^ = 0, 1, 2, 3, and keep the ground state. 
For the quench parameters used below we verified that | V'o) al- 
ways belongs to the zero magnetization sector We then diag- 
onalize the evolution Hamiltonian H using iterative Lanczos 
steps, calculating the first 500 hundred lowest energy eigen- 
states. To obtain a good approximation of the full sector we 
checked that the normalization condition was satisfied up to 
10^**, i.e. 1 — X]n=i-Pn < 10^^. The time series statistics of 
observables are obtained using 400,000 random samples in a 
time window [0, Tmax] where Tmax is chosen more than two 
orders of magnitude larger than the smallest energy gap in the 
system. 



The Hamiltonian ([T]i represents a fully interacting quantum 
system. In order to gain some initial insight into the phys- 
ical properties of this system, we numerically calculate the 
five lowest energy levels of the model as a function of the ap- 
plied local field ^,5 on the four adjacent spins (representing 
the subsystem S) in a chain of 16 spins for different ratios of 
the nearest and next-nearest neighbor couplings (see Fig.[T]i- 

Phase diagram at hs — 0: The phase diagram for the 
model at zero field and in the large N limit is well known 
and has been described for example in |22]. At zero field and 
J2 = we have a Heisenberg spin chain with only nearest- 
neighbor coupling, which is solvable using the Bethe Ansatz. 
For both positive Ji and J2 the system is frustrated. For small 
values of J2 a gapless antiferromagnetic phase is present. At 
J2/J1 « 0.241 a gap opens up, and the system remains 
gapped for all finite J2/J1, but the gap closes in the limit of 
J2/J1 -^ 00, where the model is approximately described by 
two weakly coupled chains. At zero field and J2/J1 = 1/2 
- the Majumdar-Ghosh point - the ground state of the sys- 
tem factorizes and can be determined analytically. Namely, 
the ground state manifold is two-fold degenerate spanned by 
\4'i)j \4'2) where \(j>i) consists of singlets between neighbor- 
ing dimers, and |(/)2) is |(/)i) translated by one site. 

The model in a finite local field: For small but finite local 
fields - roughly hs < 0.3 - we encounter a regime, in which 
the inter spin coupling dominates, and the local field acts as 
a perturbation. The dependence of the eigenenergies on the 
local field here is relatively flat. Degeneracies which occur 
for zero field are lifted. 

For intermediate local fields - 0.5 < hs < 2 - the Heisen- 
berg coupling and the local magnetic field compete. In this 
regime we observe numerous level crossings. 

For large fields - hs > 2.5 - the energy levels are dom- 
inated by the applied local magnetic field, and thus simply 
decrease linearly with its amplitude. This is caused by the 
alignment of the affected spins along the field direction, i. e. 
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Figure 2: Equilibration statistics for A'^ — 16, Ns ~ 4, J2/J1 = 0. 
The quench is from from hg / J\ — 3.5 to h^ / Ji — 3. In panel a) 
the probabilities spectrum p„ = p (En) is shown as a function of the 
energy. Panels b), c), d) show the full statistics of Pc, Pdsj Pms 
respectively. Both the distribution of the LE and the magnetization 
are narrow Gaussian. Moreover Pd^ indicates a small mean Ds- All 
these facts indicate standard equilibration of Gaussian type. Quali- 
tatively similar results are observed for different values of the ratio 
J2/J1. 



7715 —^ 1/2. thus maximizing the contribution of the Zeeman 
energy term, Ezeeman = —hs ■ Ns/2. In this regime one can 
treat the Heisenberg interaction as a perturbation on the field 
Hamiltonian for the spins with an applied field. In this regime, 
the slope of all the energy levels in Fig. [T] approach what one 
would expect from a dominating Zeeman term -^ — > Ns/2. 

Having these regimes of the model in mind, we consider 
the following quench scenarios: small quenches within each 
of the regimes and large quenches across different regimes. 

The simplest case is a small quench within the regime of 
large dominating local fields, where the energy levels in good 
approximation simply vary linearly with the local field ampli- 
tude hs- As an example of this case, we consider the quench 



from an initial field h 



(») 



3.5 to an evolution field h 



(/) 



for the three different couplings J2/ Ji = 0, 0.5, 1 (see Fig.|2|i. 
In all of these cases the Loschmidt Echo as well as the local 
magnetization show Gaussian distributions with very small 
variances, indicating a good, straightforward. This is exactly 
the behavior expected for small quenches in regular regimes 
[9]. The distribution of the quantity Dg in these cases also 
resembles a Gaussian, only experiencing a sHght asymmetry 
due to the non-linearity of the norm. So in regular regimes 
our example system shows good equilibration, both locally 
and globally. 

Another set of relatively simple quenches are those from 
large fields to zero field. As the system in this case is strongly 
perturbed, one would expect numerous excitations across a 
wide range of energies. According to £ = l/dP^ (p) a large 
number of excitations causes a small Loschmidt Echo. More 
specifically such quenches lead to an exponential distribution 
of the Loschmidt Echo with an average very close to zero 
Jgl]. We numerically observe such a behavior when quench- 



Figure 3: Equilibration statistics for A'' = 16, Ns ~ 4, J2/J1 — 0. 
The quench is from from hg /Ji = 3 to hj /Ji — 0. Panels as in 
figure|2] As can be seen from the log-linear plot in the inset of panel 
b), the distribution of the LE is exponential with a very small mean. 
Qualitatively similar results are observed for different values of the 
ratio J2/J1. 



ing from large to zero local magnetic fields for all three cou- 
plings. Figure [3] shows a representative example of a quench 
from hg = 3 to hj = with only nearest-neighbor cou- 
pling. In this quench, one also obtains single peaked and rela- 
tively naiTow distributions of the local magnetization and the 
quantity Dg. This indicates measure concentration and ex- 
actly the kind of strong local equilibration that is possible even 
for closed quantum systems which has been discussed in fl. 
Results for different couplings are not shown here, but very 
similar 

For small quenches in the regime of dominant Heisenberg 
couplings, we observe different types of equilibration and a 
strong dependence on the inter-spin coupling. This is dis- 
cussed in the following section. 



IV. NUMERICAL RESULTS FOR NON-TRIVIAL 
QUENCHES 

Here we investigate small quenches of the form 



ffo(4^) 



H{hi^'^) 



with h 



(*) 



0.2 and h 



(/) 



on 



systems of size iV = 16 and different coupling ratio. In 
particular we present three simulation for quenches on a 
subsystem S of four spins {Ns = 4) and coupling ra- 
tio J2/J1 = 1, 0, 1/2, and three analogous simulations 
for Ns — 3. For all the simulations the initial ground 



state is located in the 5*^^ 
magnetization. 



sector of vanishing total 



A. Quenches on four adjacent spins 

Quench 1: The first example we consider is a quench on a 
system of equal nearest and next-nearest neighbor couplings 
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Figure 4: Quench 1. Equilibration statistics for A*' 
J2/J1 = 1. The quench is from from hg / J-_ 
0. Panels as in figure|2] 
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0.2to/iyVJi 
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( Ji = J2). The system is quenched from Ho{hg = 0.2) to 

H(hg — 0). The resulting equilibration statistics are shown 
in Fig.m As one would expect for a small quench, the energy 
probability distribution p(£'„) in Fig.|4^) is dominated by the 
ground state (p{Eo) = 0.86). A more interesting feature is an 
additional sizable contribution of the first excited state. Note 
that the population of the first excited state is about two or- 
ders of magnitude larger than the population of all the others 

(p{Ei)<$:p{Ei),i>i). 

The existence of these two dominating modes leads to a 
double-peaked distribution function of the Loschmidt Echo 
Pc, which is clearly observed in Fig. |4j5). In this case, one 
can effectively neglect the weights p„ with n > 2 in Eq. |6l 
and treat the model as an effective two-state system. Hence we 
obtain the following approximation for the Loschmidt echo 



m 



pI 



2 
Pi 



-2popiCos[{Ei - Eo)t\. 



(7) 



^£ 



The corresponding probability distribution function can be 
calculated analytically: 



Pc{x)=C 



I- C{X2- Xi) 
TT^(X - Xy){x2 - X) 



(8) 



where x\ — L — 2popi, X2 = C + 2popi denote the lower 
and upper edges of the distribution. Equation ([8]l, with po and 
pi obtained numerically, is shown as a red line in Fig. ^) 
and shows excellent agreement with the full calculation. This 
result shows that, after the quench, the system oscillates be- 
tween the two lowest states of H, |0) and |1), clearly indi- 
cating absence of equilibration. Such a lack of equilibration 
can also be observed in the variance, A£ = {x2 — xi) /\/8 
which is of the same order of the support of the distribution. 
Especially since the evolution Hamiltonian is translationally 
invariant, one would also expect to see this non-equilibrium 
behavior locally. Accordingly, the statistics of Ds in Fig.HJ;) 



shows lack of equilibration. Its distribution function also has 
a large spread and two maxima, one of which shows a similar 
divergence as the Loschmidt Echo. The asymmetry and the 
lack of a second peak are not numerical artifacts but arise due 
to the high non-linearity of the trace distance, namely its ab- 
solute value. A simplified example of this behavior is given in 
the Appendix. 

It is interesting to note that the distribution of the local 
magnetization is Gaussian displaying good equilibration prop- 
erties. The reason for this becomes clear when looking at 
Eq. (|2]i. The largest (in modulus) contribution to that sum 
is given in principle by the term with n — 1 and m ~ 0, but 
this contribution vanishes as (l|O|0) = 0. Hence the sum in 
Eq. (|2]) contains a large number of terms of the same order of 
magnitude. In this situation (O (t)) can be seen as a sum of 
many independent variable converging to a Gaussian thanks 
to a central limit theorem argument. 

This situations makes clear that even when Ds is large one 
can still find observable (those with very small or zero weights 
in the low energy states) which show good equilibration prop- 
erties [i24tl . 
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Figure 5: Quench 2. EquiUbration statistics for TV = 16, Ns ~ 4, 
J2/J1 — 0. The quench is from from /ig /Ji — 0.2 to hg / Ji — 
0. Panels as in figure[2] 



Quench 2: The second quench we show is from a field of 
hg — 0.2 on four adjacent spins to zero field, but using only 
nearest-neighbor coupling, i.e. J2 = 0. The energy proba- 
bility distribution of this quench is shown in Fig.|5^). In this 
distribution the ground state is by far dominating {p{Eq) — 
0.99), the probability of the next highest populated state is 
two orders of magnitudes smaller Overall, there still is a con- 
centration of excitations in the low energies (E < —5.5). 

As indicated by the distribution of p{En), the Loschmidt 
Echo mean of this quench is much closer to one and its vari- 
ance is about an order of magnitude smaller compared with 
the case J2 = Ji (Fig. E}?)). Its distribution still shows two 
maxima, but they are less pronounced and smoother 

In general, a better approximation to the LE distribu- 
tion function can be obtained by retaining the N„iax largest 
weights Wij = piPj in equation (|6]l. This procedure will be 



performed for all of the following quenches. For the quench 
discussed here, it is enough to keep N^ax — 5. The resulting 
approximate distribution is shown as a red line in FiglSjj). 

Both the LE and Ds indicate a much better equilibration 
compared to Quench 1 . This is in contrast to the behavior of 
the local magnetization in Fig.HlJ), which is again Gaussian, 
but shows a larger rather than a smaller variance. 
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Figure 6: Quench 3. Equilibration statistics for A*' — 16, Ns = 4, 
Jij J\ = 1/2, Majumdar-Ghosh point. The quench is from from 
hf/Ji = 0.2 to h'-P/Ji = 0. Panels as in figure[2] 



Quench 3: The last quench we want to discuss in detail 
is the special case of J2 = 1/2. The evolution system in 
zero field is at the Majumdar-Ghosh point. The system has 
a two-fold degeneracy at its minimum energy, consisting of 
two states, each of them a product state of nearest-neighbor 
singlets. This degeneracy is lifted by the pre-quench applied 
local field, with hg = 0.2 acting on four adjacent spins. As a 
consequence the system cannot be treated as non-degenerate. 
To take the degeneracies into account the corresponding for- 
mulas have to be modified, e.g. p includes off-diagonal terms, 
where £"„ = £",„ for n ^ m. Furthermore, because of degen- 
eracies, the proof for inequality (|5]l given by [4] does not hold 
in this case. 

Inspecting p{En) in Fig. [6^), as in the case of J2 = 
the population of the ground energy level is by far domi- 
nating, p{Eo) = tr[IlE„Po] = 0.96 (U-Eo is the projector 
onto the bi-dimensional ground state manifold). We also ob- 
serve many more excitations, none of which is considerably 
larger than the others. The distribution is concentrated at en- 
ergies below E ^ —5, but even for E > —5 there are some 
weak excitations. Quite a few of the lowest energy eigen- 
states, are not populated at all, but protected by symmetry, 
e.g. p{Ei) ^ p{E2) ^ 0. 

Since the p(£'„) are relatively widely distributed, and be- 
sides p{Eo ) there is no other dominant contribution, we expect 
and numerically obtain a Gaussian distributed LE, as shown in 
FiglSl As in the previous quench we can approximate the LE 
truncating Eq.|6]to the first N„iax largest (in modulus) terms. 
As explained in [10], in order to recover a Gaussian distribu- 
tion we have take N^ax large, here N^ax = 20 is needed. 



The corresponding probability distribution is the line shown 
inFig.|6]3). 

Compared to the quench at J2 = the LE mean is smaller 
and the distribution shows a higher variance. So far, this in- 
dicates relatively straightforward equilibration. Therefore the 
distribution of Ds shown in Fig. |6l;) is quite surprising. It 
shows a slightly asymmetric Gaussian with the lowest vari- 
ance of all the quenches discussed, but with a relatively large 
meaning = 0.419. 

A large Ds does not mean bad equilibration for all local 
observables, but it indicates that there is at least one local ob- 
servable, which equilibrates badly. As in the other quenches, 
the distribution function of the local magnetization is a Gaus- 
sian centered around zero. Here, in agreement with the LE, its 
variance is larger than in the case J2 = 1. 



B. Quenches on three adjacent spins 
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Figure 7: Quench 4. Equilibration statistics for TV — 16, Ns ~ 3, 
J2/J1 ~ 1. The quench is from from hg /Ji — 0.2 to hg / Ji — 
0. Panels as in figure[2] 

We also show the equilibration statics of the same three 
quenches just discussed, with the same quench amplitude but 
with a magnetic field acting on three instead that on four ad- 
jacent spins, i.e. A^ = 16 and Ns = 3 (Figs. |2]|9]l. These 
quenches do not need to be discussed in the same detail, but 
indicate that some of the observed patterns are not restricted 
to Ns = 4 and provide a few interesting other features. Note 
also that the corresponding variances of the considered ob- 
servables C, Ds and ms in the case of Ns = 3 are much 
smaller than for Ns = 4. 

Quench 4: In the case of J2/J1 = 1 (Fig.|2]l one obtains a 
similar dominance of two lowest-energy modes as for the case 
Ns ~ 4. Contributions of other states are only one order of 
magnitude smaller than the two dominating modes. The LE 
distribution accordingly is double peaked but smoother and 
with a wider spread compared to the 4-site local quench that 
is about an order of magnitude smaller To obtain a good ef- 
fective description using for the LE we use here Nmax = 5. 
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plays similar features as in the case J2I J\ = 1, with all the 
distributions that tend to be broader 

Quench 6: The last quench shown in Fig. |9] is the case 
•h/ Ji — l/2> where the evolution Hamiltonian is at the 
Majumdar-Ghosh point. As in the case for Ns — 4 one ob- 
serves a relatively large number of non-vanishing modes and 
a Gaussian distribution of both the LE and the local magneti- 
zation. The approximation shown in Fig.|9}5) is obtained using 
Nmax = 10. Although the distribution of Dg is highly peaked 
with a small variance its average is relatively large, indicating 
that there can be local observables experiencing poor equili- 
bration. 



Figure 8: Quench 5. Equilibration statistics for TV = 16, Ng ~ 3, 
J2/J1 = 0. The quench is from from hg / J\ = 0.2 to hg /Ji — 
0. Panels as in figure|2] 
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Figure 9: Quench 6. Equilibration statistics for A'^ — 16, Ns = 3, 
J2/J1 — 1/2, Majumdar-Ghosh point. The quench is from from 
hf/Ji = 0.2 to h\I^/Ji = 0. Panels as in figure[2] 



The contribution of these few other states is even more visi- 
ble in the distribution of Ds, which again shows two maxima, 
but is much less spiked. Differently from the case Ns = 4, 
the local magnetization is double peaked following a distribu- 
tion function similar to the one of the LE. Correspondingly the 
variance is roughly 4 times as large as for Ns = 4, indicating 
worse equilibration. This can simply be explained by the fact 
that the two dominant modes, namely the ground state and the 
second excited state in this case have a finite cross term in the 
local magnetization. 

Quench 5: The quench at J2/J1 = shown in Fig.[8]dis- 



V. CONCLUSIONS 

In this paper we have numerically studied different local 
quantum quenches on frustrated spin chains. More precisely 
we considered both global and local equilibration indicators, 
the Loschmidt echo and the subsystem magnetization. To gain 
a better understanding of the equilibration properties of the 
subsystem we also studied the distance of the reduced density 
matrix from its time average: Ds {t) — \\ps (t) ^ As||- Our 
numerical simulations show that a variety of possible scenar- 
ios take place depending on the interplay between quench am- 
plitudes and coupling constants. In particular, for large and 
intermediate quenches the long time probability distribution 
of the Loschmidt echo resemble an exponential and a Gaus- 
sian respectively, as it was first observed in [9J- When the 
Loschmidt echo mean is very small the subsystem does equi- 
librate even when the bounds in [4] are too loose to be ap- 
plicable. On the contrary, when the Loschmidt echo mean is 
large and thus no general conclusion can be drawn a priori, 
we observe all possible situations. Namely we observe both 
cases where the subsystem equilibrates {quenches 2, 4) and 
where it does not {quenches 1, 3, 5, 6). Moreover even when 
the subsystem does not equilibrate as a whole, we still en- 
counter instances where one observable (the magnetization) 
shows good equilibration properties {quenches 1, 3). 

In conclusion, we would like to note that the identification 
and characterization of the time-scales associated with the dif- 
ferent equilibration phenomena we observed lies outside the 
scope of the long-time statistics techniques adopted in this 
paper Addressing this problem in its generality is perhaps 
the most exciting next challenge in the quest for understand- 
ing equilibration dynamics of unitarily evolving quantum sys- 
tems. 
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Appendix 

Let us assume that the system is divided into a subsystem S 
and an environment _E. The two corresponding bases are given 
by {\a) , |6) , • • • } and {\a) , |/3) , • • • }. To model quench 1 
and to simplify the calculation we further assume that only 



two eigenstates of the closed system S (E) E contribute to the 
full density matrix. (A reasonable simplification of the p„ dis- 
tribution in Fig. EH).) To model the strongly coupled system 
and to obtain a non-unitary evolution of the subsystem, we 
have to introduce entanglement. This can be done by consid- 
ering the two states 

|1) = (|a,a) + |6,/3))/V2, (9) 

|2) = i\a,a)-\b,/3))/V2. (10) 

Using the initial state |\['o) = ci |1) + C2 |2) we compute 
Ds (t) and its distribution. Given the energy difference cu = 
E2 — El of the two eigenstates we obtain 



cicl\l){2\e' 



h.c. 



(11) 



Tracing out the degrees of freedom of the bath this simplifies 
to 

Ps-ps ^ Trs (p - p) (12) 

= cic*{\a){a\-\b){b\)e''^' + h.c.. (13) 



Identifying C1C2 = y/piP2e^'^ and using a matrix representa- 
tion we get 



PS - PS == \/PlP2 cos {ujt + if) 



We finally obtain 



Dsit) = y/PlP2 \C0S {ujt + ip)\ 



(14) 



(15) 



Its distribution function can be calculated analytically, giving 
a simple description of a divergence similar to what we ob- 
serve numerically, 



Pds{x) 



lim i f 5{Ds{t)-x)dt (16) 



= — / ^ WpiP2 |cos {ujt)\ - x) dt (17) 



2/n 



\JP\P2 - X'^ 



(18) 



Note that in this case Ds takes values in Ds G [O, ^^plpi] , 
so such a P [ds) has only one square root singularity at the 
upper edge. If we simply take the two largest amplitude con- 
tributions of Fig. [4^), we can estimate the peak in Fig.|4};) to 
be around \/0.86 • 0.13 = 0.33. Given the great simplification 
this estimate is quite useful. 



